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We construct finite mass, asymptotically flat black hole solutions in d = 4 Einstein- 
Q \ Yang- Mills theory augmented with higher order curvature terms of the gauge field. They 

possess non-Abelian hair in addition to Coulomb electric charge, and, below some non-zero 
^Jy critical temperature, they are thermodynamically preferred over the Reissner-Nordstrom 

solution. Our results indicate the existence of hairy non-Abelian black holes which are 
stable under linear, spherically symmetric perturbations. 

> 

Introduction.— In recent years it has been realized that the electrically charged Reissner- 
Nordstrom (RN) black hole, when considered as solution of a more general theory, may become 
unstable to forming hair at low temperatures. This has lead to the discovery of some holographic 
models for condensed matter systems, and, in particular, to a gravitational description of 
superconductivity (see [1] for a review). 

The case of Einstein- Yang-Mills (EYM) model with negative cosmological constant A in 
d — 4 spacetime dimensions provides an interesting illustration of these aspects. As shown in 
[2], there is a second order phase transition between the RN-anti-de Sitter solutions, which 
are preferred at high temperatures, and symmetry breaking non-Abelian black holes, which are 
preferred at low temperatures. In [2], A plays an essential role; although electrically charged 
hairy black holes do exist also in a Minkowski spacetime background |3J, they have rather 
different properties as compared to the anti-de Sitter (AdS) solutions in [2]. In particular they 
do not emerge as perturbation of the RN black holes, and, similar to the well-known d = 4 
asymptotically flat, purely magnetic EYM solutions [I], are also perturbatively unstable. 

However, one might take the view that in the strong coupling regime the (A = 0) EYM 
theory is incomplete. Perhaps the simplest possibility to describe this situation is to supplement 
the action of the EYM model with higher order curvature terms, for both gravitational and 
gauge field sectors. As discussed e.g. in the inclusion of (string theory inspired-) corrections 
to the gravity action does not lead to qualitatively new features. By contrast, we will here argue 
that the situation is different vis a vis the inclusion of higher order Yang-Mills (YM) curvature 
terms. This possibility has been overlooked so far in the literature. The first relevant order in 
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this case is the fourth, in which case the most general such density added to the Lagrangian 
consists of the four terms, 

C s = Cl Tr {F^F pa F^F^} + c 2 Tr {F, u F^F pa F^} (1) 
+ c 3 Tr {F^F^F^F^} + c 4 Tr { F ilv F vp F p \F Xil } , 

with some constant coefficients Cj. A particularly priviledged such combination, which we 
adopted here, is that with c\ = c 2 = — 4c3, C4 = 0. In that case, C s features only the second 
power of any "velocity field" and is a causal density just like the Gauss-Bonnet term in gravity 
[6] or the Skyrme [7] term of the 0(4) sigma model. With this specific choice of the constants 
Cj, the Lagrange density (JT|) is nothing else than the trace of the square of the 4— form cur- 
vature F pvpa = {F^ v , F pa ]}. This is the second member of the YM hierarchy [8], providing 
a natural generalization of the usual YM model. A convenient way to express this system is 

Tr /(F^F^) 2 j, where a tilde denotes the Hodge dual. 

Notwithstanding our specific choice for the constants c, in (III) , we have verified that for 
certain other choices, some salient features of the solutions discussed in this work, in particular 
the instability of the RN black hole, persist. 

The model.— Ignoring for simplicity other possible corrections, we consider the following 
action for the model 

S = J dV=s[~^ - ^Tr {F, v Fn + ^Tr {(F^) 2 } ] , (2) 

(here we have set 47rG/e 2 = 1, such that the only parameter of the theory is r). 

In what follows, we shall prove that the presence of the last term in (j2J) leads to an instability 
of the RN black hole, together with the occurance of stable black holes with non-Abelian hair 
outside the horizon. We shall restrict attention to the following spherically symmetric Ansatz: 

dv 2 

d s 2 = ^+ r 2 {d6 2 + sin 2 9d<f) 2 ) - Na 2 dt 2 , (3) 

where N, a are functions of r and t in general. The minimal gauge group for which the super- 
position of a Coulomb field and a non-Abelian hair is not forbidden by the 'baldness' theorems 
[9] is £77(3). Then, as in the r = case in [3], we shall restrict to an £77(2) x £7(1) truncation 
of the £77(3) group, the general spherically symmetric ansatz for the gauge potential being 

A = { (uT 3 + UT 8 )dr + (wT t + wT 2 )d9 + ((wT 2 - wT x ) sin 9 + cos 9T 3 ) d<\> + (vT 3 + VT 8 )dt} , (4) 

where u, w, w, v and U, V are functions of (r, t) and Tj are the standard generators of the £77 (3) 
Lie algebra. 

For static solutions, one can set the functions u, w, v and U to zero without any loss of 
generality, resulting in the equations 

, Ar /2 (1-w 2 ) 2 r 2 V' 2 (l-w 2 ) 2 V 2 , 2a , 2 

m> = Nw' 2 + + — — + t± ^ , a' = —w' 2 , 

2r la r A a z r 

„ f N' a\ , w{l-w 2 ) 2r{w 2 -l)V' 2 n 
Jy a r z i\ r z Na z 
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together with the first integral for the electric potential, 



a_( 2r(l - 



v ' = Q— ( 1 + —2 ) • (0) 

with Q an arbitrary constant. 
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The RN solution corresponds to a vanishing SU(2) field, w(r) = ±1, and m{r) = M — 2 

2 
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er = 1, V(r) = § — Q/r. This solution has an outer event horizon at = M + a/ M 2 — Q 
which becomes extremal for Q = M. 

Solutions with nonzero magnetic gauge fields should also exist. However, one can see that, 
for Q 7^ 0, the first integral © excludes the existence of particle-like configurations with a 
regular origin. Thus, the only physically interesting solutions of this model describe black 
holes, with an event horizon at r = > 0, located at the largest root of N(rh) = 0. The 
regularity conditions at the horizon imply the following series expansion there 

m \ T ) = TT + m n r ~ r h) + ■ ■ ■ > a \ r ) = °h H \r -r h ) + ..., 

2 r h 

w(r) =w h + wi(r - r h ) + . . . , V(r) = v x (r - r h ) + . . . , (7) 

wtiere v x - r 4 +2r(1 _^ )2 , mi - ' Wl ~ — (i-2 mi )r h ^ — ■ i1] 1S also 

straightforward to show that the requirement of finite energy implies the following asymptotic 
behavior at large r 

O 2 J 2 JO 

m {r) = M —-£- + a{r) = 1 - — + ... , w{r) = ±1 + - + ... , V{r) = $ - ^ + . . . . (8) 

Once the parameters an, Wh an d J, M, Q are specified, all other coefficients in ([7j), (jHJ) can 
be computed order by order. M and Q correspond to the mass and electric charge of the 
solutions; other quantities of interest are the Hawking temperature T H = j-cr(rh)N'(rh), the 
entropy S = = i:r\ and the chemical potential $. J here is an order parameter describing 
the deviation from the Abelian solution. 

The results.— The solutions of the equations (JSJ), ([6]) interpolating between the asymptotics 
©, ([H]) were constructed numerically, by employing a shooting strategy. Finite mass, non- 
Abelian black holes exist for any r > 0. For given Q, rh, t, the solutions are found for discrete 
values of the parameter Wh, labeled by the number of nodes, n, of the magnetic YM potential 
w(r). To simplify the picture, in this work we have restricted attention to solutions with a 
monotonic behavior of the magnetic gauge potential w(r) (and thus n = 0, 1 only). 

In characterizing the non-Abelian configurations, it is convenient to introduce the quantity 

« = (9) 

For < k < 1/2, the solutions can be thought of as nonlinear superpositions of the RN and the 
purely magnetic SU(2) black holes in [10J. In particular, they are unstable since the magnetic 
gauge potential has always n > 1 nodes. Moreover, their free energy F = M — T^S is always 
greater than that of the RN configuration with the same T# and Q. 
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Figure 1: The profiles of a typical non-Abelian nodeless solution are shown together with the 
corresponding potentials for the perturbation equations ( IT2l) . (Tl4"l) below. 



The picture is very different for k> 1/2. In this case, for given k, we notice the existence 
of a set of non-Abelian solutions emerging as perturbations of the RN black holes, for a critical 
value of the charge to mass ratio. This instability of the Abelian configuration is found within 
the Ansatz for values of the magnetic gauge potential w{r) close to the vacuum values ±1 
everywhere, w(r) = ±1 + eW(r). The perturbation W(r) starts from some nonzero value at 
the horizon and vanishes at infinity, being a solution of the linear equation 

<w)'-<i-^=o, do, 

where N = (1 — — The second term in this equation can be seen as an effective 

mass term /i 2 for W near the horizon, with /x 2 ~ 1 — 2«;(Q/r/ i ) 4 . We have found that for any 
k > 1/2 (i.e. fi 2 < 0), an instability occurs for a critical value of the mass to charge ratio of 
the RN solution. An approximate value of this ratio is found by using an asymptotic matching 
expansion for the approximate solutions of (TTOT) at the horizon and at infinity, the result 



Q _ 2V6V1 + Vl + 48^ 
M 7+ VI + 48k 

providing good agreement with the numerical data. 

This unstability signals the emergence of a symmetry breaking branch of the non-Abelian 
solution bifurcating from the RN black holes. In contrast to the solutions with k < 1/2, here 
we notice the existence of a fundamental branch of solutions without nodes in the magnetic 
gauge function w(r), see e.g. the typical profile shown in Fig. 1. 

Among the solutions bifurcating from RN black holes, those with 1/2 < k < 2 are of special 
interest, sharing some striking similarity with the picture found for the EYM-AdS system in 
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[2] . A plot of the horizon function of the temperature reveals the existence of a single 

branch of non-Abelian solutions, which exist below a critical temperature only (the range 
of the scaled temperature tjj — T^Q being < tjj < 0.0203). In a canonical ensemble, the non- 
Abelian black holes within this range of k exist for a finite interval of (i.e. of the entropy) only. 
These solutions always possess a positive specific heat, the Hawking temperature vanishing for a 
minimal value rj i " Mn ' ) = \p2r\l yplk — 1 of the event horizon radius. As r — > r\^ m \ an extremal 
non-Abelian black hole solution with a regular horizon is approached, the charge to mass ratio 
of this configurations being always greater than one. Furthermore, it turns out that the free 
energy of a RN solution is always larger than the free energy of a non-Abelian solution with the 
same temperature and electric charge. Therefore, for 1/2 < k < 2 these non-Abelian black holes 
are preferred. The difference of free energies scales like (T^ — Th) 2 near the transition point, 

signaling that this is a second order phase transition, while J ~ \J 1 — Th / T^ . Moreover, for 
the same values of the mass and electric charge, the RN solution has a smaller event horizon 
radius (and thus a smaller entropy), than the non-Abelian black hole. 

The picture is somehow different for k > 2. Again, one finds a single branch of solutions 
emerging as a perturbation of the RN black hole, nodeless configurations existing also in this 
case. However, these solutions exist for values of the temperature greater than the critical value 
Tjj and are thermally unstable. 

Some of these features are shown in Fig. 2 where we plot the scaled horizon area an = 
A H /Q 2 as a function of the scaled temperature tn = T H Q for several values of the ratio t/Q 2 . 
The branch of RN solutions is also shown there. In Fig. 3 the scaled free energy f = F/Q 
is plotted as a function of the scaled temperature tn for several values of k. The inlet there 
shows the behavior of the parameter J which enters the asymptotic behavior of the magnetic 
potential, as a function of the ratio Th/T^. 
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Figure 3: The scaled free energy / = F/Q is plotted vs. the scaled temperature tu = TuQ 
for several values of k = r/Q 2 . The inlet shows the behavior of the order parameter J which 
enters the asymptotics of the magnetic potential. 



Existence of stable solutions.— An outstanding question now is whether the F 4 term 
leads also to stable non-Abelian black holes. The fact that, for a range of k, we have found 
nodeless solutions which are thermodynamically favored over the RN black holes suggests a 
positive answer to this question. 

For simplicity, we consider linear, spherically symmetric perturbations only. Even in this 
case, the analysis is highly involved and the details will be presented elsewhere. Here we briefly 
outline just the key features. 

As usual, all field variables are written as the sum of the static equilibrium solution whose 
stability we are investigating and a time dependent perturbation. Choosing a gauge such 
that U — v — 0, one finds that the fluctuations decouple in two groups. Sw(r,t), SV(r,t), 
5a(r,t) and 5m(r,t) form even-parity perturbations, whereas 5w(r,t) and 5u(r,t) form odd- 
parity perturbations. The linearized equations imply that 8a(r,t), 8V(r,t) and 8m(r,t) are 
determined by 8w(r,t) = wi(r)e~ lQt , leading to a single Schrodinger equation 

cPw 

1 + U even (p)w 1 = tt 2 w 1 , (12) 



dp 2 

(where dr /dp = Na) with a potential 
Na 



U e ven 



(, 2rV ' 2 \ fo 2 , 8W, 2 1N Aw' 2 (l-w 2 ) 2 
1 j- 3w 2 - 1 + (w 2 - 1) K — '- 



(13) 



r2 ,' 2(1 -w 2 ) 2 r 2 V' 2 \ 32t 2 w 2 (1-w 2 ) 2 V' 2 
" ■"' { r 2 + a 2 ) ' r.-r ; (l • ^(1 - ,r 2 )^ 



which is a regular function in the entire range — oo < p < oo. The corresponding analysis 
for the odd sector is much more evolved. After much algebra, the perturbation equations for 
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Sw(r,t) = Wi(r)e tnt and 5u(r,t) = vi(r)e int can be cast in the form 

- + U odd (P*)V = A 2 *, (14) 

where ^ = [(l + 6r(l-w 2 ) 2 /r 4 )^-12r(«; 2 -l)w' /r^w^Fu dr/dp* = Na/(l + 12rNw' 2 /(r 2 (l + 
6r(l - w 2 ) 2 /r 4 ))) 1/2 and i>i(r) = F 2 \I> + F 3 ^f'. The potential U od d and the functions Fj have 
complicated expressions depending on the equilibrium functions m, a, w and V^, with F\ > 0. 
For nodeless solutions, U Q dd and F are regular in the entire range — oo < p* < oo. 

For both equations (112j) and ( fT4"l) . the potential vanishes near the horizon and at infinity. 
Then standard results [TT] imply that there are no negative eigenvalues for Q 2 (and hence no 
unstable modes) if the potentials U even and U dd are everywhere positive. 

Our results indicate that this is indeed the case for some of the solutions with k > 1/2, see 
e.g. the inlet in Fig. 1. Interestingly, approaching the extremality appears to imply generically 
positive values for the potentials in f|T2|) . f }T4"|) . 

Therefore we conclude that, in contrast to all other known d = 4 asymptotically flat hairy 
black holes with non-Abelian gauge fields only [4] , at least some of our solutions here are linearly 
stable. 

Further remarks.— We close with some remarks on the generality of the results in this 
work. First, we remark that the F 2 and F 4 terms in ([2]) are the pieces which also enter the 
Lagrangian of the non-Abelian Born-Infeld theory [T2] describing the low energy dynamics of 
D— branes. Although the equations of motion of that model coupled with gravity are more 
complicated and do not admit the RN black hole as a solution, it is natural to expect that the 
picture we have found here will share some similarities with the results in that case. 

Non-Abelian fields featuring both the F 2 and F 4 terms appear also in the higher loop cor- 
rections to the action of the d — 10 heterotic string [13]. However, a generalization in that 
framework of the solutions considered here is not an easy task, due to the occurrence there of 
a variety of other fields. Previous work in this direction [5] indicates that the solutions with a 
standard F 2 term only, share the basic features of the EYM black holes in [10], in particular 
being unstable. It appears however, that inclusion of higher order gauge field curvature terms 
could lead to a very different situation (for example, we have found rather similar results to 
those discussed above, for a generalization of ([2]) with an extra dilaton field, a Gauss-Bonnet 
term and gauge group 5*0(5)). Thus we expect the following picture to be generic: the F 4 
term introduces a supplementary interaction between electric and magnetic fields, which, for 
some range of the parameters, implies a tachyonic mass for the vacuum perturbations of the 
non-Abelian magnetic fields around the Abelian solutions. Then the Abelian gauge symmetry 
is spontaneously broken near a black hole horizon for some critical value of the charge to mass 
ratio, with the appearance of a condensate of magnetic non-Abelian gauge fields there. The 
possible relevance of these aspects in providing analogies to phenomena observed in condensed 
matter physics is yet to be explored. 
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